High fidelity state transfer in binary tree spin networks 
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Quantum state propagation over binary tree configurations is studied in the context of quantum 
spin networks. For binary tree of order two a simple protocol is presented which allows to achieve 
arbitrary high transfer fidelity. It does not require fine tuning of local fields and two-nodes coupling 
of the intermediate spins. Instead it assumes simple local operations on the intended receiving node: 
their role is to brake the transverse symmetry of the network that induces an effective refocusing 
of the propagating signals. Some ideas on how to scale up these effect to binary tree of arbitrary 
order are discussed. 

PACS numbers: 03.67.Hk, 03.67.Lx 



I. INTRODUCTION 

The paradigmatic approach to quantum communica- 
tion assumes the possibility of "loading" quantum infor- 
mation (i.e. qubits) into mobile physical systems which 
are then transmitted from the sender of the messages 
to their intended receiver. Such flying quhit architec- 
ture for quantum communication has found its natural 
implementation in optics where photons play the role 
of information carriers. In many respect this appears 
to be the most reasonable choice, specially when long 
distance are involved in the communication. However 
the recent development of controllable quantum many- 
body systems such as optical lattices [l|, phonons in ion 
traps [2] J Josephson arrays i^, and polaritons in opti- 
cal cavities (J], makes it plausible to consider alternative 
quantum communication scenarios such as the so called 
quantum wire architectures [^ . Here the transfer of quan- 
tum information proceeds over an extended network of 
coupled quantum systems (e.g. spins) which are at rest 
with respect to the communicating parties. In this case 
the messages are encoded into the internal states of the 
spins while the information flow proceeds by their mu- 
tual interactions which, when properly tuned, induce a 
net transfer of messages from two separate regions of the 
network [a, |7|, la] . The quantum wire architecture is of 
course of limited application, since it assumes the sender 
and the receiver to have access to the same quantum 
network (in any real implementation the latter will al- 
ways have a reduced size). However these techniques 
may play an important role in the creation of clusters of 
otherwise independent quantum computational devices. 
Furthermore the study of quantum network communica- 
tion protocols is an ideal playground to test and device 
new quantum communication protocols. 

Perfect transfer among any two regions of a quantum 
network can always be achieved if one allows the commu- 
nicating parties to have direct access on the individual 
nodes of the network (for instance this can be done by 
swapping sequentially the information from one node to 
subsequent one) . These strategies are however extremely 
demanding in terms of control and, even in the absence 
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FIG. 1: Spin tree of second order. The nodes are connected 
through edges which describe the XY interactions defined 
by the Hamiltonian ^, and are identified by a double index 
(a, h), with (0, 0) corresponding to the leftmost graph element. 
Inset: Auxiliary spin added to the BT in order to homogenise 
the effective couplings in the block form representation ((6)1. 
The additional spin provide also a controllable trigger to start 
the information transfer. 



of external noise, are arguably prone to error due to the 
large number of quantum gates that have to be applied 
to the system. A less demanding approach consists in 
fixing the interaction of the network once for all and let- 
ting the Hamiltonian evolution of the system to convey 
the sender message to the receiver. In this context per- 
fect transmission can be achieved either by engineering 
the spins couplings [3, [13, [ill, [l4| , or by choosing proper 
encoding and decoding protocols [Tj, [Tj, [H, [la, [l3| • 

In this paper we discuss the propagation of quantum 
information over Bifurcation Tree (BT) quantum net- 
works. Together with the star configuration the BT 
configuration is arguably the most significant network 
topology in circuit design. The former are typically used 
as hubs to wire different computational devices (for an 
analysis of such system in the context of spin network 
communication see Ref. [11]). Star configurations have 
been also extensively studied for entanglement distribu- 
tion [131 and cloning [19]. BT networks instead are em- 
ployed to route toward external memory elements (i.e. 
database). The information flow on unmodulated and 



uncontrolled BT was first discussed in Ref. [20| while, 
more recently, BT quantum networks have been em- 
ployed to design efficient quantum Random Access Mem- 
ory elements (2ll|. 

The paper is organised as follows. In Sec. [IT] we start 
analysing the first non trivial BT system introducing the 
notation and setting the problem. In Sec. IIIII we then 
describe a transfer protocol that allows one to deliver a 
generic quantum message to any desired final edges of the 
second order BT network by exploiting simple end gates 
operations. In Sec. lIVI we discuss various techniques that 
allow us to scale up the protocol adapting it to BT of 
arbitrary order. The paper finally ends with the conclu- 
sions and discussion in Sec. Ivl 



II. SYSTEM DESCRIPTION 

First order BT networks are just particular instances 
of star networks [ll|. Consequently the simplest non- 
trivial examples of BT networks is the second order one 
shown in Fig. [1] In the following we will assume the lines 
connecting the nodes to represent XY (exchange) spin 
interactions (the results however can be generalised to 
include XX Z or Heisenberg couplings). The resulting 
Hamiltonian is thus 



H 






E^(^ 



1) 



(1) 



where the summation is performed over all couples i,j 
which are connected through an edge, where (^l^y^^ ^^^ 
the Pauli matrices associated with the i-th node, and 
where the w^'s appear in consequence of the interaction 
with local magnetic fields (in this expression the label 
i an j stands for the joint indexes (a, 6) of Fig. [1]). As 
usual [5| we assume that initially the system is in the fer- 
romagnetic "all spin down" ground state |0) = |0 • • • 0). 
At time t = we then place an (unknown) qubit state 
\tp) = a|0) +/3|1) on the left-most site (0, 0) (for instance 
by swapping it from an external memory). With this 
choice the global state of the network is now described 
by the vector 



l*i 



a|0) 4-/311(0,0) 



(2) 



where |l(a,f))) represents the network state where the node 
(a, b) is in the spin up state |1) while the remaining ones 
are in the down state |0), i.e. 



|l(a..)>^|0- 



•0 1(,,,)0- 



•0) 



(3) 



Knowing that the z component of the total spin is pre- 
served by the Hamiltonian evolution of the system (i.e. 
[H, 5*^°'] = 0) we can conclude that the dynamics is cos- 
trained in the subspace of single-flips: On this subspace 
H acts in a very simple way that can be inferred from 
the graphical structure of the network, i.e. 
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(a,b)) 



■'{a,b)\^{a,b) 
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|l(c,d)) , 



(4) 



where the sum is taken over all sites (c, d) connected with 
(a, 5). Our goal is to find a procedure that would allow 
us to transfer the qubit state jV') to one the right-most 
sites (2, b) with b € {1, 2, 3, 4} of our choice, i.e. 



I*i 



I* 



(b)\ 
fin/ 



a|0)+/3|l(2,b)) 



(5) 



Following Refs. fT(j. ^T] one could try to solve this prob- 
lem by fine tuning the parameters Jo and loj of H in such 
a way that the free Hamiltonian evolution of the system 
will be able to transform I^Pin) into |^gn) after some time 
interval r j22l |. This however is in general a quite com- 
plex calculation which entails to solve an inverse eigen- 
value problem. Moreover, if any, the solutions obtained 
using such strategy will be arguably highly asymmetrical 
in the distribution of the local magnetic fields (jJj's. To 
avoid all this, here we will pursue a different approach 
by limiting the freedom one has in choosing the Hamil- 
tonian parameters but, as in Refs. [IJ, [la)ll3> allowing 
local manipulation on the receiving node of the network 
(i.e. (2,6)). Under these conditions we can show that a 
simple protocol exists that realises the transformation ([3]) 
with arbitrary accuracy. It assumes an homogeneous net- 
work structure where all the ratios loj/ Jq are chosen to be 
identical and equal to some fix value, and it is composed 
by the following three-steps: 

1. the system is allowed to evolve freely under the 
action of H for some time r; 

2. at this point on the receiving node (2,6) is per- 
formed a fast (ideally instantaneous) local trans- 
formation PS(2.b); 

3. the network is then let evolve for an extra time 
interval 2r. 

During the first step, due to the homogeneity of the 
Hamiltonian, the information flows along the left-right 
axis of the network while delocalising along the south- 
north axis. The value of r is approximatively the time 
interval an excitation takes to travel from the leftmost 
node (0,0) to the rightmost column formed by (2.1), 
(2,2), (2,3) and (2,4). The role of the local transfor- 
mation PS(2,6) of step two is to brake the south-north 
symmetry of the resulting state by flipping the sign of 
a specific wave vector component. The system is then 
let evolve freely for a time interval which twice the ini- 
tial one: this is approximatively the time it takes an 
excitation to leave the rightmost column, " bounce back" 
the leftmost node, and return to the rightmost network 
column. Due to the symmetry brake introduced at the 
second step, however, the signal will now not diffuse over 
all the four sites (2.1), (2,2), (2,3) and (2,4), but in- 
stead it will focus on the intended receiving node (2,6). 
A detailed description of the protocol will be presented 
in Sec. Hm 



A. Diagonalisation of the Hamiltonian 

To solve our problem we can exploit the fact that the 
ground state |0) of the network does not evolve to re- 
strict ourselves to the case a = 0, i.e. |V') = |1)- We then 
simplify the structure of the Hamiltonian (1) assuming 
all ujjS to be identical, i.e. luj = uj. In dealing with mag- 
netic spins this means we are applying an homogeneous 
magnetic field of constant strength B oc — a; all over the 
system. We could set uj — 0, since the energy is defined 
up to a constant, but we let it be nonzero to guaran- 
tee that |0) is the ground state of the system. We now 
choose the following basis for the single excitation sector, 
that divides the Hamiltonian in invariant blocks: 
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f4)^ 43(11(1,1)) -|1(1,2))) 



\V5) = 5(|1(2,1)) + |1(2,2)) - |1(2.3)) - |1(2,4))) 

)) 



\vi) 
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H2,l) 



(|1(2,1)) 



II 



K2,2); 

(2,2)) 



II 



K2,3); 

(2,3)) 



L(2,4) 



|1(2,4))) 



In this basis the matrix representing H is given by 



\ 




LlI J 
J UJ 



UJ 
UJ 



(6) 



/ 



with J = \f2jQ. This shows that the evolution of the 
network can be effectively described as three indepen- 
dent linear chains, the first composed by 4 nodes, and 
the other by 2 elements each. The basic idea in deriving 
the above basis is that any state in the form |l(a,6)) is de- 
coupled from the "singlet" superposition ^(|10) — |01)) 
of the two nearest neighbours qubits on its right. E.g. 
The two states 4R(|l(2,b)) — |l(2.b+i))) with fe = 1,3 are 
decoupled from the whole network and provide an alter- 
native basis for the block ^3. 

Of special interest for us is of course the block Bi which is 
the only one to have an overlap with the input state ([2]). 
It is clear that the case Jq = J would be much simpler 
to deal with (in this case for instance one could adapt 
the linear chain analysis of Refs. [y, [l^l to simplify the 
calculation). Such option however is not possible if we 
assume the coupling strengths of the network to be fixed 
a priori. Anyway we can use a trick to obtain the same 
result without adjusting the coupling strength which, as 
discussed in the final paragraph of the present section, 
allows us to improve also the controllability of the setup. 
Suppose in fact to modify the system by adding an ad- 
ditional spin connected only to site (0, 1) with the usual 



XY coupling of strength Jq, as shown in the inset of 
Fig. [TJ With this choice the Hamiltonian ([1]) is replaced 
by 

//new ^ i/ + t «ux<l + ^fux<l) + f «ux + 1)- (7) 

Now that we enlarged the Hilbert space we have to deal 
with the 9-dimensional space of single flips. However 
since the singlet state -75(11(0,0)) ~ |laux)) is decoupled 
from the rest, if we encode the "logic" state |1) on the 
sending end of the network as |wo°™) = -75 (|l(o,o)) + 
|laux)) instead of using |l(o,o))j not only we recover a 
dynamics costrained in an 8-dimensional space, but we 
obtain also an effective coupling of strength J = v^ Jo 
between |w5°*) ^^^^ I'^^i) (here |laux) is the analogous of 
the states ^ with the spin up localised on the auxil- 
iary node). The four dimensional block of our effective 
Hamiltonian thus becomes: 



^(4)^ 




(8) 



Following Ref. [lC>| this can be easily put in diagonal form 
obtaining the eigenvalues 



-fci = 2 



E- 



_ 2uj+{V5-l)J 



E4 



■2 — 2 ■ 

_ 2uj+{V5+l)J 



(9) 



with the corresponding eigenstates described by the vec- 
tors 



|ei 



_1 l + \/5 _l + \/5 1 
' 2 ' 2 ' 



e2)^^i=(l,i-^ i-^ 

63 



_ 1 



1, 



V5-V5 



_-, 1-V5 1-V5 1 1 
' 2 ' 2 ' / ' 



|e4) 



ysWi 



^ i+VE i+VE ^ 



(10) 



(expressed in the basis |wo°*), |«i), |«2) and l^s)). 

A part from simplifying the spectral properties of the 
Hamiltonian, the introduction of the site aux adds also an 
additional feature that substantially enhances our abil- 
ity of controlling the system. We have already men- 
tioned the singlet state -7s(|l(o,o)) ~ |laux)) is decou- 
pled from all others vectors of the system. Therefore 
we can "entrap" our qubit of information at the left- 
most end of the network for as much time as we like 
by encoding its logic 1 component in such singlet state. 
When we want the transfer to start we simply apply a 
local Phase shift PSaux on the auxiliary spin that in- 
duces the mapping |laux) — ^ — |laux)- This will transform 
^(|l(o,o))-|laux)) into Iv^""") bringing the encoded mes- 
sage into the four-dimensional subspace associated with 
the Hamiltonian ([5]) and allowing the first step of the 
above protocol to begin. 



III. TRANSFER PROTOCOL 

In this section we analyse in detail the performance of 
the protocol defined in Sec.|TTl Without loss of generality 
we consider the case in which the receiving node is (2, 1), 
i.e. 6=1. We recall that our aim is to obtain the 
transition |t)JJ™) 



|l(2,i))i Ei-iid we notice that 



|l(2,l)) = 5(l^3> + |w5> + |f6> + |«7» 



(11) 



The protocol: step 1. In the first stage of the proto- 
col the system is initialised into \v{)) and freely evolves 
for some time t. The goal here is to find uj and r such 
that this vector is mapped into jua), which represents a 
symmetric combination in which the input excitation is 
spread all over the rightmost nodes of the network. As 
already noticed, this process is formally equivalent to the 
information transfer along a linear chain of 4 spins cou- 
pled by uniform XY first neighbours interactions. From 
Ref. [10| we know that such transferring cannot be exact. 
Nevertheless the transfer fidelity can be made arbitrar- 
ily close to one. Indeed defining U{t) = exp[— iiJr] and 
using Eq. Q we have 



(^3|f/(r)|«r")=5]e-^'^*(«3|e,)(e.|i;i) 



(12) 



k=l 

%/5 — 5 { „~iEiT ^ — iEiT\ _|_ \/5+5 ('^— i-E2T 

20 



-(e 



) + ^IP(e 



-iE'iT 



This will be exactly one if one could find r such that 

g-iiSir ^ g-iSar ^ _i g^^^ ^^iE2r ^ ^-lE^r ^ ^ ^^^^ 

though these conditions are impossible to be satisfied ex- 
actly .lOy an approximate solution is obtained by choos- 



ing 



and 



7+x/5 7 _ 7+\/5 7 



2r7, + l 



J 



-TT , 



(13) 



(14) 



with n integer. Under this assumption Eq. (J12p yields 

(z;3|C/(r„)|i;S™) = i(l+e-''^"), (15) 



where (^„ = y/b{2n + l)7r. The exponential in Eq. p^ 
never takes the value 1 but, since v5 is an irrational 
number, approaches it indefinitely. Therefore for any e > 
we can choose n such that |(w3|[^(T„)|wg°^)| > 1 — e. 
As a result the state of the system, with high accuracy, 
is now described by the vector {v^). 

The protocol: step 2. As a second step we act lo- 
cally on the node (2, 1), applying the phase shift unitary 
transformation PS21 which changes the sign to the state 



1 



(2,1)/ 



PS21II 



211-1(2,1)/ 



II 



(2,1)/ 



(16) 



while preserving the remaining single excitation states. 
This can be done, for instance, by acting with an intense 
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FIG. 2; Minimum "infidelity" {1 ~ F„) achievable according 
Eq. (O for n G [0, N] as a function of iV, i.e. (1 - F)min = 
niin„^jv(l — Fn). Notice that already for n = 8 we get values 
of F greater than 1 — 10~®. From the plot we can infer an al- 
most linear dependence of (1 — F)niin with respect to A'', yield- 
ing the following (approximated) behaviour Fmin ~ 1 — cN ' 
with 7 ~ 1 and c being constant. 



magnetic field which acts locally on (2,1) for a time in- 
terval shorter than the characteristic times of the Hamil- 
tonian H. When acting on jws) the unitary PS21 yields 
the following transformation. 



PS2ib) = i(|f3)-k) 



\V6) 



\V7)) 



(17) 



This superposition contains the four states that compose 
the state |l(2,i)), but the relative phases are wrong - see 
Eq. (|lip . Luckily third step fixes this issue with free 
evolution only. 

The protocol: step 3. Finally we just have to wait for 
a time 2t„ and the relative phases adjust themselves to 
give the state |l(2,i))- In fact by explicit calculations it 
can be shown that 

(i;3|C/(2T„)|t;3> = i(l + e-^2^")-l, 

(18) 
{vk\U{2Tr,)\vk) = -e-^^- ~ -1 , for fc = 5,6,7. 

Both expressions are justified by the fact that e~"^" ~ 1, 
and they imply that after third step we have reached 
state |l(2.i)) with as good approximation as we like. 

The above operations can be summarised by the ap- 
plication of the unitary operator 

K = t/(2T„) PS(2,1) f/(T„) . (19) 

Therefore the resulting transfer fidelity can be expressed 
as 

_.-,. 2 

,,ncw\ 1 2 



F„ 



1(1(2.1) iKl^-o"™) I' 



1-He- 



3e-''^" -f 



8 



(1 



-tifrx 



= i cos2(^„/2)[3 - cos2((^„/2)]2 , (20) 

1 giving us the desired 



which approaches 1 if e *'''" 
result - see Fig. [5J 
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FIG. 3: Modified tree for the purpose of scaling up. Each 
rightmost qubit is accompanied now by an auxiliary, moreover 
all the rightmost couplings are adjusted by a factor l/\/2. 



IV. SOME IDEAS TO SCALE UP THE SYSTEM 

Unfortunately our protocol doesn't extend easily to 
higher-order trees (or at least we couldn't find a sim- 
ple way of doing it). The idea we pursued in trying to 
scale up a second order BT is to connect in some way the 
ends of a tree to the beginning of another. The result- 
ing structure isn't anymore a tree of the type described 
above, but it's still a valid mean to obtain a larger num- 
ber of outputs. As an example we could connect (say) 
two second order trees to the ends of a first order tree 
to obtain a 8-outputs quantum switch, or four second 
order trees to the ends of a fifth second order tree to 
obtain a 16-outputs quantum switch - see Fig. [3] The 
former setup can be solved by properly merging the pro- 
tocol of Ref. |ll| with our second order BT propagation 
scheme: this however will require to employ non-uniform 
magnetic fields at least for the first spins and does not 
admit simple concatenation. We thus decided to focus 
on the second architecture which instead can be trivially 
concatenated to form larger setup. We found a relatively 
simple way to make the required connections, but at the 
expenses of considering some coupling strength engineer- 
ing and including antiferromagnetic interactions, which 
means that the "all down" configuration is no longer the 
ground state, although still stationary. A combination 
of time evolution and a phase shift will do all the work. 
First of all each receiving end of a tree must be accompa- 
nied by an auxiliary qubit, as done before for the sending 
end. In analogy to what happened introducing (0, 0)aux, 
it is easy to see that the rightmost singlets 



\S,,) ^ ^(|1(2,6)) - 11(2,=^.))) , 

are isolated, while the corresponding triplets 

|i26)^^(|l(2,.)> + |l(2"r^))), 



interact with the network and evolve, with an effective 
coupling strength -\/2 times the original one. In order 
for our protocol to be still valid, we need to modify the 
coupling strengths of the rightmost branches so that ma- 
trix ([8]) (now with Jo = J) remains unchanged. Moreover 
the local operation on the receiving end (2, b) nrust now 
be performed simultaneously on (2, b) and (2, 6)aux, i-e. 
we must now apply PS26 -1- PS2fc''. In this way once the 
excitation reaches one of the end-triplets of the tree it 
can be trapped there with a Phase Shift on the auxil- 
iary qubit of that site, storing information in the relative 
singlet. To clarify this, we outline that in this new con- 
figuration our protocol is capable of achieving the (ap- 
proximate) transfer 



«|0>+/3|l(o,o)> 



Now by applying a local phase shift PS2b 
transformed into 

a\0)+f3\s2b), 



a\0)+p\t2b). (21) 

''"'' the state is 



(22) 



which is decoupled from the rest. 

If by some means we could transfer this state to the sin- 
glet at the beginning of the next tree (denoted by primed 
indexes), i.e. obtain the state 

a|0)-f/3i=(|l(o.,o')>-|l(o'>)))^"|0)+/'l^o'O'> , (23) 

we could then perform the local operation PSfp'^^o') to 
obtain the corresponding triplet state 

a\0)+l3\to,o'), 

that can be transferred along the new tree with the usual 
protocol. The structure shown in Fig. 0] (that we will call 
"singlet link" ) achieves perfect transfer between two sin- 
glets, since in the subspace {|s2r), |lo), |so'0')} il is equiv- 
alent to a chain of length 3 with constant couplings ^1C|, 
moreover the evolution of |i2r) is decoupled from that of 
|io'0') thanks to the opposite signs of the couplings along 
the vertical axis. The lines stand for XY interaction of 
strength J (ferromagnetic) and —J (antiferromagnetic) 
respectively. As an example we have considered site (2,1) 
of a second-order tree plus its auxiliary qubit, connected 
with site (0',0') plus its auxiliary of another tree. We 
outline again that we are working in the subspace of sin- 
gle flips, as our Hamiltonian still conserves 5*^°*. In the 
considered example we have 

H\s2i)=u;\s2i) + V2J\lo), 
i7|i2r)=co|t2i) + J|l(i,i)), 
H\sq>q>) =w|so'0')+V2J|1o), 

i7|io'o'> =w|io'o'>+%/2|l(oM')); 

where |io'0')^^(|l(o',o')> + |l(o'>)))- 
We can see from the above equations that once the in- 
formation enters a tree through a triplet state it doesn't 




FIG. 4: The "singlet link" used to connect the outputs of the 
first tree to the inputs of the next column of trees (only the 
first one is represented here). The auxiliary qubits are not 
represented for simplicity. 



come out of it until we make a Phase Shift on the desired 
end (and at the right time!). At this point the informa- 
tion goes to the singlet and propagates to the starting 
singlet of another tree, thanks to the singlet link, then 
it is transferred to the corresponding triplet with a local 



Phase Shift and propagation begins on the next tree. We 
shall repeat this procedure until information reaches the 
desired end on the last array of trees. Of course we must 
control a priori the total error due to the presence of sec- 
ond order trees, so to fix a "single-tree time" f which 
gives a satisfactory overall transfer fidelity. 



V. CONCLUSIONS 

In this paper we have presented a protocol for quan- 
tum state transfer on BT spin networks of order two. As 
in Ref. [ll, [la [i3 it is based on the local operations 
which must be performed on the receiving nodes. Dif- 
ferently from |16l . llTJ I however it does not involve swap- 
ping operation between the receiving nodes and external 
memories and arbitrarily high fidelity can be obtained 
in just three operational steps. Generalisation of this 
techniques to higher orders BT is currently under inves- 
tigation: arguably this will involve more complex ends 
gates operations possibly on more than one of the right- 
most nodes. We have however provided a simple way to 
scale up the problem by concatenating smaller BT net- 
work through connecting gates which can be turned on 
and off by simple local phase gate transformations. 

We thank R. Fazio and D. Burgarth for comments and 
discussions. This work was in part founded by the Quan- 
tum Information program of Centro Ennio De Giorgi of 
the Scuola Normale Superiore. 
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